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A derivation of the t — J model of a highly-correlated solid 
is given starting from the general many-electron Hamiltonian 
with account of the non- orthogonality of atomic wave func- 
tions. Asymmetry of the Hubbard subbands (i.e. of "elec- 
tron" and "hole" cases) for a nearly half-filled bare band is 
demonstrated. The non-orthogonality corrections are shown 
to lead to occurrence of indirect antiferromagnetic exchange 
interaction even in the limit of the infinite on-site Coulomb 
repulsion. Consequences of this treatment for the magnetism 
formation in narrow energy bands are discussed. Peculiarities 
of the case of "frustrated" lattices, which contain triangles of 
nearest neighbors, are considered. 



The problem of strong correlations and magnetism in 
many-electron (ME) systems is one of the most impor- 
tant in the solid state theory. Derivation of the simplest 
ME model describing these phenomena in the case of 
narrow energy bands (strong Coulomb interaction) was 
proposed by Hubbard [Q. This model was studied in a 
great number of papers. 

A detailed investigation of fcrromagnetism in the Hub- 
bard model in the limit of the infinite on-site Coulomb 
repulsion U was performed by Nagaoka He proved 
rigorously that the ground state for the simple cubic and 
bcc lattices in the nearest-neighbor approximation with 
the number of electrons A^e = A'' ± 1, N being the num- 
ber of lattice sites, possesses maximum total spin, i.e. 
is saturated ferromagnetic, since this ordering provides 
the maximum kinetic energy gain for an excess electron 
(hole). (The same statement holds for the fee and hep 
lattices with the transfer integral t < 0, — N + 1, 
or t > 0,Ne = N — 1.) The picture of saturated fcr- 
romagnetism is preserved at small, but finite concentra- 
tions of current carriers c = \Ne/N — 1|. In the case of 
a half-filled band (A^e = N), \t\ -^U the ground state is 
antiferromagnetic because of the Anderson's kinetic ex- 
change (superexchange) interaction ^ which occurs in 
the second order in |t|/C/. This interaction is due to the 
gain in the kinetic energy at virtual transitions of an 
electron to a neighbor site, which are possible provided 
that the electron at that site has an opposite spin direc- 
tion. In systems with a large finite U and ^ N, a 
competition between ferro- and antiferromagnetic order- 
ing occurs. Nagaoka ||] has put forward a criterion of 
fcrromagnetism, which was based on the condition of the 
spin-wave spectrum stability (the magnon spectrum can 
be calculated also from the spin Green's function, see. 



e.g., Ref. [Q). This criterion has the form 

ac > \t\/U 



(1) 



where the constant a of the order of unity depends on the 
lattice structure, a — 0.246 for the simple cubic lattice. 
At the same time, pure antiferromagnetism is stable at 
Ne = N only, and a phase separation ||] takes place 
provided that the condition (|^) is violated. 

A version of the Hubbard model in the limit U ^ oo 
with inclusion of the antiferromagnetic exchange inter- 
action J, the t — J model, is widely applied now to de- 
scribe copper-oxide high- Tc superconductors and related 
systems (see the review [^). The quantity J is usually 
related to the Anderson's exchange (J = — 2<^/C/), but 
sometimes it is considered as an independent phenomeno- 
logical parameter. Recently the difference between the 
t~ J model and large-?/ expansion of the Hubbard model 
has been discussed 0. 

In his pioneering work jl] Hubbard neglected the differ- 
ence between the atomic functions and orthogonal Wan- 
nier functions and did not treat the non-orthogonality 
problem. At the same time, in the limit of large U it is 
more natural to use well-localized atomic wavefunctions. 
In the present paper we develop a consistent formulation 
of a general ME model with strong correlations with the 
use of the Bogoliubov approach 1^ and discuss the elec- 
tron spectrum picture and magnetism formation in this 
model. 

We start from the general Hamiltonian of the ME sys- 
tem in a crystal 



H 



2 |r,-r,- 



where 



Ze2 



= 5]i;(r-R,) 



(2) 



(3) 



is the periodic crystal potential. For simplicity, we con- 
sider a non-degenerate s-type band (the case of degener- 
ate bands, originating from the atomic d- or /-functions, 
can be considered in a similar way ||,|l§). To pass to 
the second quantization representation we have to use 
orthogonal wave functions. However, the atomic wave 
functions 



(4) 



(s is the spin coordinate) do not satisfy this condition 
at different sites v. We apply the orthogonalization pro- 
cedure developed by Bogoliubov H, which enables one 
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to construct the perturbation theory in the overlap of 
atomic wavefunctions belonging to different lattice sites. 
To lowest order in this parameter the orthogonalized 
functions read 

^•^ (r) = (r) - ^7 ^ ip^, (r) (5) 

where the sum goes over the nearest neighbors, 

l = lu^'^ / drLpl,{r)ip^{r) (6) 



is the non-orthogonality integral. Substituting (|^) into 
and using the orthogonality of the spin wave func- 
tions Xcris) we obtain the ME Hamiltonian of the "polar 
model" 



(7) 



Vi V2 1^3 ^4 '^tl CTl '^Viai ^1^4 <^2 '^f 3 



Then we derive 

W = e ^ (fT, a) + ^ X, (2, 2) 

C'^' 0)^-2 (-^, 2) + X,, (2, ~a)X,, (0, a)]} 

)]} (12) 

where 

rti, ^^^n^a = al^CLya ^ ^^[Xy{a, a) + X,^(2, 2)], 



and 



U — 1 1, 1/1/1/ ■) Q 



(13) 



(14) 



are the Hubbard parameter and the Coulomb integral at 
different lattice sites, 



where 



(8) 



is the one-electron level in the central potential of a given 
site w(r). 



rfrCi(r)(-^A + y(r))V,,(r) (9) 



J dY(j)l^ {v)v{r ~ Ri.^ (r) 



are the transfer matrix elements between the sites Ui and 

' drdr'Ci (r)V:2 (^')^i:^^-3 (r)V'i/. (r') 

(10) 

are the matrix elements of interelectron Coulomb repul- 
sion. We have neglected in (||) and (||) the influence of 
potentials of the sites v' ^ v and v' ^ i^i , respec- 
tively, since the corresponding terms contain extra fac- 
tors v(v — Ry/)-i/;i,(r) with v' ^ v which are small due to 
the decrease of the potential w(r) with increasing r (in 
other words, the crystal potential is small in the region 
between the lattice sites). We have also dropped in (|^) 
small integrals which contain v(v — Riy')|'0i,(r)|^. 

Further we retain in (^ only one- and two-site terms 
and pass to the representation of the Hubbards operators 
Xv{y^ , A) ll^ which transform the state | A) (A = 0, cr, 2) 
at the site v into |A'), 



a|:„ =X,((T,0)+aX,(2,-f7) 



(11) 



^(00) _ ^ 

''1/11/2 ~ ''1/11/2 

+ (22) _ . + 

''1/11/2 — ''1/11/2 ' 1/1 1/1 1/2 1/1 

,(02) _ ,(20) _ , , r 

''1/11/2 ~ ''1/11/2 ~ ''1/11/2 ' 1/1 1/1 1/2 1/1 



(15) 
(16) 
(17) 



are the transfer integrals for empty states (holes) |0) and 
doubly occupied states (doubles) |2), and the integral of 
the double-hole pair creation. 

Estimations of the parameters of the Hamiltonian (|^) 
with non-orthogonality being neglected are presented in 
Ref. jl]. In the expressions for e, f/ and Q the non- 
orthogonality corrections are small in the overlap and 
do not play a role. On the other hand, at calculating 
other parameters in (|l^) we have to take into account 
the second term in (^). We obtain for the integral of the 
"direct" exchange 



d If 11/2 -^1/1 1/21/21/1 



(18) 



JviV2 271/^2^2^1/11/2 + ^{U + Ql/il/2 )7l/il/2 



where 



Jl'll/2 -^1/1 1/2 1/2 1/1 ' -^1/11/2 — -^1/1 1/1 1/2 1/1 7 



(19) 



the integrals / being calculated for the atomic functions 
(f. The expressions (|l5|)-(|l7|) take the form 



+(00) _ + 
1/1 1/2 ~ ''1/1 1/2 



^Li/l ^1/11/2 ~^ 2Z/,yj,y2 {U + Ql/ll/2)')'l 



1/1 1/2 — ''I/II/2 ^ ^ 



V1V2 2 (^71/11/2 + ' 



(20) 
(21) 

1/1 1/2 )7i/ii/2 (22) 



All the terms in (|8|), as well as in (^)-(^), are of the 
same (first) order in the overlap. Moreover, the quan- 
tity U"f should be larger than J/7 and \t\. Indeed, 
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in the case of narrow bands the interelectron repulsion 
(which determines U,Q, J and L) and the crystal poten- 
tial (which determines e and t) are of the same order of 
magnitude, although the crystal potential is expected to 
be somewhat larger. At the same time, as follows from 
@i (P^, t (or Q, J and L) contain another small factor 
that was considered above: these integrals include the 
product of the function (j)v{r) by the potential v{r — R^/) 
(or the electron repulsion) corresponding to other lattice 
site, which "cuts" a rather small region r ~ R^/. Thus 
the "on-site Coulomb" (non-orthogonality) contributions 
should dominate in ( ^l|) and (|2^). 

It should be noted that in fact the transfer integrals 
( pl| ) and (|2^ ) are to be calculated with the use of many- 
electron wave functions (see the review ^) which are, 
generally speaking, not reduced to the Slater determi- 
nants and not factorized into one-electron ones. For ex- 
ample, the general Hartree-Fock approximation in the 
atom theory (see Ref . ) uses the radial "one-electron" 
wave functions which depend explicitly on the ME atomic 
term jP). The transfer integrals are expressed through 
the corresponding ME wave functions as M 



,(rr',r"r"')= J l[d{r.,s,}%^r^t 

f V(r, 



1^2 r" 



(23) 



2m 



Therefore the integrals (P0D-(p2[) can be different even 
at neglecting interatomic Coulomb interactions and non- 
orthogonality. 



The electron spectrum of the model (12) in the sim- 
plest "Hubbard-I" approximation (which corresponds 
to a "mean-field approximation" in the electron hopping, 
the on-site Coulomb repulsion being taken into account 
in the zero-order approximation) is given by B 



+ (24) 



T^{[tr'M- t 



f)(l-(n_.))-C/]^ 



+4(tf))"(n.)(l-(n_.))}V2 

One can see that, unlike the standard consideration [Q, 
the Hubbard subbands (^^ turn out to have quite dif- 
ferent widths even in the paramagnetic case for a nearly 
half-filled band ((n+) = (n_) ~ 1/2). In particular, in 
the large-C/ limit we have 



and non-orthogonality integral (||) . Thus an appreciable 
asymmetry of the cases Ne < N and > N can occur, 
the bandwidth in the case of hole conductivity being con- 
siderably smaller than in the electron (double) case. This 
circumstance may be important, e.g., for copper-oxide 
high- Tc superconductors. Of course, more advanced ap- 
proximations than ( p^ ) should be used in the antiferro- 
magnetic state, especially in the two-dimensional case. 
In the latter situation, the low-energy electron spectrum 
is determined by the scale J rather than t (see Ref. [p^). 

Generally speaking, we have to take into account in 
( |T2| ) also three-site "operator" Coulomb contributions to 
the transfer integrals. 



The quantity I^^ ^ 



(26) 



1/1 1/' 1/21/' is small in comparison with 
LviU2 — IviuiW2Ui due to the decrease of the Coulomb 
interaction with distance. A peculiar situation occurs in 
the case of a "frustrated" lattice where equilateral trian- 
gles of nearest neighbors are present (e.g., the triangular, 
fee and hep lattices), so that the site v' can be the near- 
est neighbor for both the sites vi and V2. Then substitut- 
ing ( [l0| ) into (|2^) yields the non-orthogonality correction 
of order of C/7. Such corrections are important for the 
calculation of electron spectrum (e.g., in the "Hubbard- 
1" approximation) and yield a suppression of the above- 
considered asymmetry of the hole and double subbands. 
In such situations, additional three-site "exchange" terms 
also arise , so that the non-orthogonality problem p^ ] 
is more difficult. 

The kinetic exchange (Anderson's superexchange) in- 
teraction occurs in the second order in t'^u^X ■ Performing 
the canonical transformation which excludes the double- 
hole pair creation and annihilation terms from the Hamil- 
tonian ( p^ ) we derive 



= ^ 



(' 



(02) 



U 



-{4(S,,S,, 



(27) 



As follows from (g2|) , the nominator in (gTj) is determined 
not only by the bare hopping, but also by the Coulomb 
interaction. Formally, the kinetic exchange interaction 
survives even in the limit U ^ 00 owin g to the non- 
orthogonality contributions. Combining ( p7|) and ( [T^ ) 
we obtain the expression for the total effective exchange 
parameter in the case of a half-filled band (c = 0) 



(1) 

kcr 

(2) 



(1 



T7^\^) _ _ I -1(22, 



- {n-a))t 

M + u 



(00) 
k ' 



(25) 



so that, according to (20),(|2l|), the bare hopping in- 
tegral t determines the bandwidth of "holes", and the 
bandwidth of "doubles" is mainly (in the above-discussed 
sense) determined by the intrasite Coulomb interaction 



Jeff = J -2 it 



'-(02) 



= J--i'^Q + 27< - 2(< + Lf/U 



(28) 



Note that the terms of the order of ^"^U are canceled 
in Jeff- First three terms in ( |2^ ) coincide with the 
corresponding result for the two-site problem (hydrogen 
molecule) jl^] and yield an antiferromagnetic exchange 
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interaction. As mentioned above, the crystal potential 
should be somewhat larger than the Coulomb interac- 
tion, and in the large-i7 limit the main contribution to 
Jeff reads 



eff 



-2j\t\ 



(29) 



where t is assumed to be negative. We see that in the case 
Ne < N the ratio of Jeff to bandwidth is proportional 
to the overlap parameter (rather than to \t\/U <C 7 as 
in the standard consideration). High values of the Neel 
temperature, which are typical for the layered copper- 
oxide systems, may be related to this fact. 

By analogy with the consideration of Ref. ||] we have 
the criterion of ferromagnetism 



2ac|i(^^)| > -J, 



eff 



(30) 

> N. Under 



where A = for TYg < iV and A = 2 for Ne 
the assumption 7 ^ |t|/C/ ^ 7^ the criterion (30) takes 
the form 



ac > 



-Jeff/m) 

-Jeff/hU) 



7, Ne<N 

2\t\/U, Ne>N 



(31) 



which is quite different from (0) for Ne < N. In the case 
Ne > N the result (|3^) is formally similar to (|l|) , but has 
a different origin. 

We see that a strong indirect antiferromagnetic inter- 
action occurs in narrow energy bands, the ferromagnetic 
exchange owing to the motion of current carriers being 
more strongly suppressed in the case of "hole" conduc- 
tivity {Ne < N) than in the "electron" case. The situa- 
tion may change considerably in the case of degenerate 
energy bands where the indirect interaction can be ferro- 
magnetic owing to the intraatomic Hund exchange ||l9|,^ . 
Note that many-band effects are often assumed to be es- 
sential also for the usual itinerant magnetism (see, e.g., 
Ref. [0). ^ 

From the experimental point of view, the narrow-band 
ferromagnetism is not a too wide-spread phenomenon. 
It takes place, e.g., in the systems Fei_a; 002,82 0, 
C0S2, Cr02 However, the degeneracy of conduction 
band plays an important role in the electron structure 
of these systems. At the same time, ferromagnetism is 
not observed in the copper-oxide systems. The above- 
discussed modifications (in comparison with the original 
Hubbard's treatment ||]) in formulation of the simplest 
s-band model and similar considerations of more realistic 
models (see Refs. []6|,pO[) may be useful for explaining this 
fact. 
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